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Abstract. 

We present a proposal for realizing local decoherence-free evolution of given 
entangled states of two two-level (TL) ions. For two TL ions coupled to a single heavily 
damped cavity, we can use engineering reservoir scheme to obtain a decoherence- 
free subspace which can be nonadiabatically controlled by the system and reservoir 
parameters. Then the local decoherence-free evolution of the entangled states are 
achieved. And we also discuss the relation between the geometric phases and the 
entanglement of the two ions under the nonadiabatic coherent evolution. 
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1. Introduction 

It is well known that one of the major obstacles for practical manipulation of quantum 
information is the effect of decoherence due to coupling to environment. So how 
to protect quantum information from these effects become a major challenge. A 
number of strategies have been proposed to fight decoherence, such as quantum error 
correction [TJ [2J E] , decoherence-free subspace (DFS) |H El El [7] , dynamical decoupling 
(DD) P, El HU], engineering reservoir [TT], etc. The general framework for DFS was 
introduced by Zanardi et al. [I] in a spin-boson model undergoing both dephasing and 
dissipation. In the presence of environment, DFS is a set of all states which are immune 
to decoherence processes induced by interaction with bath. These are groups of states 
that have robust symmetry properties. The quantum information can be encoded in 
DFS in order to increase reliability of quantum information processing [T2l [T3l E]. 
On the other hand Carvalho et al. have presented a proposal for protecting states 
of a trapped ion against decoherence, based on engineering of pointer states [IS]. By 
controlling the reservior, they have applied an indirect control on the protected states, 
and have shown how to protect squeezed states, approximate phase eigenstates, and 
superpositions of coherent states. Recently, Prado et al. have shown how to protect 
a nonstationary superposition states of a two-level (TL) ion. Working with an ion 
trapped inside a dissipative cavity they have constructed two classes of decoherence- 
free evolution by broadening engineering reservoir scheme [16]. Under the assumption 
of a squeezed engineered reservoir, Carollo et al. have proposed a way to observe the 
adiabatic geometric phase acquired by a protected state evolving coherently through 
the adiabatic manipulation of the squeeze parameters of the engineered reservoir [17J. 
Prado et al. have generalized the method proposed by Carollo et al. [17] and shown 
how to observe the geometric phases acquired by the protected nonstationary states 
even under the nonadiabatic evolution |16j . Authors in Refs. [15,16] only considered 
one particle system, then a natural question arises, how about multi-particle system? 
Can we use this engineering reservoir scheme to protect entanglement? In this paper 
we will answer these questions. 

In this paper we introduce a scheme, which generalizes the engineering reservoir 
scheme presented for a single ion in Ref. [IE], for protecting given entangled states 
of two TL ions by building a class of periodic nonadiabatic coherent evolution. We 
consider a system consisting of two TL ions coupled resonantly to a dissipative cavity 
and simultaneously driven by classical fields. For this system, we can obtain a time- 
independent master equation by the techniques of engineering reservoir. Then we can 
show that the given initial entangled state can be projected into a DFS, and in particular, 
modifications of the parameters of the reservoir may result in a controlled time evolution 
of the protected subspace as a whole. In this way the protected entangled state lying 
in this subspace evolves periodically. Finally, we calculate the geometric phase of the 
whole system and the subsystems under the nonadiabatic coherent evolution. 

The paper is organized as follows. In Sec. 2, we introduce the model and deduce the 
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time-independent master equation by engineering reservoir. Then we obtain the DFS 
and show how to protect given entangled states of two TL ions. In Sec. 3, we calculate 
the geometric phase of the whole system and the subsystems under the nonadiabatic 
evolution and give the relation between the geometric phases and the entanglement of 
the two TL ions. The conclusions are given in Sec. 4. 

2. Model and results 

We consider a system consisting of a pair of two identical TL ions A and B coupled 
resonantly to a single dissipative cavity mode of frequency uj with equal coupling strength 
g. In this paper we suppose that there are no direct interactions between the two ions. 
The ith ion is driven by two classical fields of frequencies uj\ with coupling strengths Vt\ 
(i = A, B; I = 1, 2) where V and '/' represent the ion and the classical field respectively. 
Within the rotating-wave approximation and a sufficiently small Lamb- Dicke parameter 
(to keep the motional state almost unchanged), the Hamiltonian modeling the system 
is given by 

H = Wfl ta + f {erf + of) + [(E w n?eWt-«t%? g 

(1) 

+ E*=i )2 n?eW?-«f%* + ga{a% + <)) + H.c], 

where a) (a) is the creation (annihilation) operator of the cavity mode, is the 
transition frequency of the ions, 4>f and 4>f are the dephasings of the classical fields, 
and a l kl = \k)i(l\, k and I being the ground state g and excited state e of the ions. In 
the interaction picture, the Hamiltonian can be written as: 

Hi = E^C^^'^ + nf^-^^) 

(2) 

+ H.c.] + [^- iAt a« + + H.c], 

where A = ui — loq and A \ = uo\ — u>o. The first and second terms of the Hamiltonian 
represent the classical field driving and the coupling between the ions and the cavity 
respectively. We suppose that the cavity mode is heavily damped with decay rate n and 
the ions spontaneously decay with rates 7*. Then the master equation describing the 
system is given by 

p = -i [H u p] + KD [a]p + Y, YD [<] p, (3) 

i 

where the Lindblad decoherence superoperator D[a]p = a pa) — a)ap/2 — pa)a/2 and 
D [o~ l ge ] p = ag e pa l eg — al g a l ge p/2 — pal g a l ge /2 describe cavity and ionic decays respectively. 

Next we engineer the appropriate interaction between the ionic levels and the cavity 
mode [16]. Concretely, we adjust the classical field u{ in resonance with the ionic 
transition frequency u>o, i.e., A\ = 0, and the detuning of the field u l 2 to A l 2 = —20^. 
By applying the unitary transformation R(t) = R A (t)®R B (t) to Hamiltonian (2), where 
R\t) = exp[-ifiit(e i ^a* +H.c.)]exp{-iff 2 t[cos(^)a*+isin(^)( e - i ^4 e -e i ^a* )]} with 
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ip % = (f>\ — 02> an d adjusting Q\ = — |A 2 ^> Q l 2 = —A ^> g, we can obtain the effective 
Hamiltonian in the rotating-wave approximation 

H2 = §[aVt + O + a(of + + A)], (4) 

where = cos(<//2)|e)i + ie~ i<ft sin(y/2)|#)i and = ie i ^sin(^/2)|e) i + 

cos(</?7 2 )b)i- 

Now applying the unitary transformation R(t) to the original master equation (3) 
we obtain 

p> = -i[H 2 , p'] + kD [a] p' + J2 1*D [° n ] P', (5) 

i 

where p' = R(t)^pR(t) and o' % = R(t)^ a l ge R(t) . In the limit where the cavity decay rate k 
is much larger than the other relevant frequencies, the cavity mode can be adiabatically 
eliminated and we can obtain a time-independent master equation just for the ions [H] 

pf 



> = YD[ai_ + al_]p' + £ ?D \a n ] p' , (6) 



where T = g 2 / k is the coupling strength of the engineered reservoir. Furthermore, under 
the assumption that V is much larger than the spontaneous emission rate Y of the ions, 
i.e., T 3> 7\ Eq. (6) becomes 



pi = TD[o-i_ + o B + _}p>. (7) 

From Eq. (7), the DFS is composed of all the eigenstates of the operator cr^L + cr^_ 
with zero eigenvalue and it is easy to prove that it is spanned by the following 
orthonormal states: the product state |1) = |+)a|+)b and the maximal entangled state 
|2) = ^(l— )a\+)b — \+)a\—)b)- It is easy to see that the state 

|* r > = v / T^ : |l) + v / ^e^|2), (1 > r > 0,p G (0,2tt]), (8) 

or a mixture state of |1) and |2) is the equilibrium state of the master equation (7). If 
the two ions are initially prepared in these states, they would remain the same. However 
if the ions starts outside this subspace, the situation gets more complicated. Because 
Eq. (7) is symmetric with respect to exchange of the two ions, one could analyze the 
system in terms of the antisymmetric 

, 9 v _ {\-)a\+)b-\+)a\-)b) 

and symmetric 

|1> = \+)a\+)b, 
|3> = \-)a\-)b, 

|4) = (\-)a\+)b + \+)a\-)b) 
a/2 
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subspaces [TS]. Generally the initial state can be decomposed into symmetric and 
antisymmetric components. The symmetric part of the initial state evolves toward 
the state |1) and the antisymmetric component |2) remain the same. 

In this paper, for simplicity we suppose that the protected state is pure. Using 
concurrence as measure of degree of entanglement between two qubits [20], we can 
obtain that the entanglement degree of the state |\& r ) is equal to r. Since we obtain the 
master equation (7) through the unitary transformation R(t), the protected equilibrium 
state \^ r ) is nonstationary in the original interaction picture. Reversing the unitary 
transformation R(t), the state \^f r ) (for simplicity we take Slf = Slf = OijO^ = Sl B = 
Sl 2 ), written in the interaction picture is 

|tf (t)) = R(t)\V r ) = Cl {t)\e) A \e)B + c 2 (t)\e) A \g) B 

(9) 

+ c 3 (t)\g) A \e) B + c 4 (t)\g) A \g) B , 
where Ci(t), c 2 (t), c 3 (t) and c 4 (t) are time- dependent parameters: 

Cl {t) = - ij ^f[e^ cos(^ - Q 1 t) sin(^ - fi x t) 
_ j<t>t sin (^ _ Q lt ) cos(^ - firf)] 
+ e- i2Q2t VT^Tcos(^- - Sl t t) cos(^ - Slit), 

C2 (t) = _^[cos(^ - firf) cos(^ - Slit) 

+ e i(^-0f ) sin (*£ _ Q^) sin (^ _ Q lt )j 

+ ie -i(^f+2f2 2t ) v ^r7 C0S ^ _ sin( ^ _ ^ 

<*(*) = ^[cos(^ - Slit) cos(^ - Slit) 

+ gi(*f-^) sin (^l _ ft^) sin (^ _ fi lt )] 

+ ie- [ ^ +2n ^VT^Tsm(^- - Slit) cos(^ - Slit), 
c 4 (t) = -^[e-tf sin(^ - Slit) cos(^ - Sl x t) 

- e" 1 *? cos(^ - Slit) sin(^ - fi^)] 

_ e -i(^+0f+2n 2 i) v / r ^: sin ^ _ Qit j gin ^ _ Qit y 

The f2 2 , 05., and ip l = <p\ — <f) 2 are adjustable parameters of the classical fields. The 
nonstationary protected state \^{t)) is allowed for a nonadiabatic coherent evolution, 
which can be manipulated through those parameters fix, f^, 01 an d { P l - Under the 
assumption that Sli=NSl 2 (where N is an integer and N ^> 1), the evolution is periodic 
and the period is equal to 7r/fi 2 . That is to say, we can obtain the protected initial 
state at time n-n jSl 2 (n = 1,2, •■•) and the concurrence of the system is invariable 
in the evolution. Concretely, suppose the initial state be |\I/(0)) = cx(0)|e) A \e) b + 
c 2(^)\e) A \g) b + c 3 (0)\g) A \e) b + C4(0)\g) A \g) b, where Cj(0) are known complex constants 
and |q(0)| 2 = 1 (i — 1, 2, 3, 4). Let t = in Eq. (10), we can obtain 4 equations, and 
then we can obtain the required parameters of the engineering reservoir <ft\ and (p l and 
the coefficients /1 and r of the corresponding state |^ r ) in the DFS, i.e., we can choose 
these parameters 4>\ and <p l to project the given initial state |^/(0)) to the corresponding 
state |\l/ r ) in the DFS. As a simple example, we suppose that the initial state is a maximal 
entangled state |\l/(0)) (Bell states). From Eqs. (9) and (10), we obtain the parameters 
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Table 1. parameters of the engineering reservoir and the corresponding state | 'J re- 
mittal state 1^(0)) parameters of the engineering reservoir corresponding state \^f r ) in DFS 

^{\e) A \e) B + \9)a\9)b) 7$ = f , if = f , 9 A = 0, V B = n ^(Ha\+)b \+)a\-)b) 

±(\e) A \e) B - \9)a\9)b) <f>? = 0,<f> B = 0,^ = 0,<p B = tt ^(\-)a\+)b \+)a\-)b) 

-L=(\e) A \g) B + \g) A \e) B ) cf>f = ~\ , <f>? = 0, V A = tt, V B = tt ^(\-) a \+) b _ \+) A \-) B ) 

-L(\e) A \g) B -\g) A \e) B ) <tf = -§,<£f = -§, ^ = tt, ^ b = tt G±(|-)a|+)b - |+)a|-)b) 



of the engineering reservoir and the corresponding protected state \^ r ) in the DFS (in 
table 1). Now we consider a more general example, |^(0)) = m\e)A\e)B + nel9 \g) a\q) b, 
where m, n > satisfying m 2 + n 2 = 1, and 9 G (0, 27r]. From Eqs. (9) and (10), we can 
obtain the parameters of the classical fields 

^ = f [1 + sign(n - m)] - arctan ( jbgjj) , 1 (n) 
V S = I t 1 + sign(ra - m)] + arctan (^j£^J , j 

and the corresponding state \^f r ) = ^Jmn e'a (|— }a\+}b — \+)a\—)b) + 
\/l — 2mn \+)a\ j t)b i n the DFS. It must be noted that although many methods have 
been proposed to protect the entanglement from dissipation, our approach is differ- 
ent. Most of the schemes to protect entanglement is static, but ours is dynamic. More 
specifically the protected state goes through a cyclic nonadiabatic coherent evolution, 
but the degree of entanglement does not evolves and remains the same. In this way we 
can change one entangled state into another entangled state with the same degree of 
entanglement against dissipation. 

All our discussions were based on Eq. (7), where spontaneous emission effects 
were neglected. However, spontaneous emission is the fundamental limiting factor for 
the existence of entanglement in a system of ions. The Eq. (6) describes the ionic 
system including the effect of spontaneous emission which introduces a coupling between 
the symmetric and antisymmetric subspaces. Next, we will analyze the spontaneous 
emission effects in our protected schemes. As a example, we consider that the two TL 
ions are initially prepared in 

\* E ) = ^=\g) A \g)B + \{\e) A \g) B - \g) A \e) B ), (12) 

Without spontaneous decay, from Eqs. (jUJ) and (10), we can see that by adjusting the 
parameters of the classical fields 

</ =<p B = vr, = 0, 0f = tt, (or < = tt, <f>* = 0), (13) 
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the initial state \^e) will undergo a coherent local evolution and the can be 

written as 



Now instead of Eq.(7) we numerically solve Eq.(6), which include the effect of 
spontaneous emission for the given initial state \^e)- We can compute the fidelity 
F = Tr[|^ E (t))(^ E (t)|R(t)p / E (t)R t (t)], where p' E {t) is the solution of Eq.(6). Within 
the regime fii = 10fi2=100g, g=500'j A , ^ A =^ B and n=3g, the fidelity F as a function 
of time is shown in Fig. 1. From Fig. 1, we can see that the spontaneous emissions 
move the system away from the protected state (^(t)). If the spontaneous emission 
rate is very small, the fidelity reduces slowly, for example, when t=1007r/f2 2 > the fidelity 
is around 96.9%. Therefore, as long as the spontaneous emission rate 7, is much smaller 
than the other relevant frequencies of the problem, we can neglect the spontaneous 
emission effect within a finite time, such as t <C 1007r/f22 in Fig.l. It is worth stressing 
that our protected scheme might be realized experimentally. The setup of two atoms 
equally coupled to a cavity mode with possibility of individual addressing has already 
been demonstrated in [21J. The large cooperativity parameter ( T = g 2 /n ^> 7) has 
been obtained in a variety of recent experiments [22j [23| 124"] . 

3. Relation between entanglement and geometric phase of the system 

From above discussion, we can see that the protected two-ion state (9) is 'dynamic' 
rather than 'static'. If Qi=NQ,2 (where N is an integer and N ^> 1), the system will 
undergo a cyclic coherent evolution. After a cyclic evolution, the system returns to its 
original state but may acquire a geometric phase. To compute the geometric phase we 
use the definition given in Ref. [25]. From Eq. (Q, after a cyclic evolution, r = 7r/f2 2 , 
the acquired geometric phase is 



From Eq. (15) it can be seen that, the geometric phase of the whole system is only a 
simple linear function of the entanglement degree r, and have nothing to do with other 
system parameters. Next we calculate the geometric phase of the subsystems and to 
study the relation between the geometric phase of the subsystems and the entanglement 
degree. Generally speaking, the state of the subsystem is no longer a pure one, so we 
adopt the definition of geometric phase for mixed states under bilocal unitary evolution 
[26] . If the Schmidt coefficients are nondegenerate, after a cyclic evolution, the geometric 





(14) 




(15) 




Protecting entangled states of two ions by engineering reservoir 



8 



phase of the subsystem can be written as 

P\t) = axgEiLiPkO/klt/^K) 

(16) 

xeM-r Q (Pk\U l \t)U\t)\^ k )dt)l 

where p k is the Schmidt coefficient, \fik) is the corresponding eigenstate of the reduced 
density matrix p l (obtained after tracing over the other ion) and U l (t) is a local unitary 
evolution operator acting on the i-th ion. 

For our system, if the initial state is not maximally entangled (r ^ 1), the Schmidt 
coefficients are nondegenerate. Using Eq. (16), we can obtain the geometric phase of 
the subsystems under bilocal unitary evolution R{t) 

(3 A = (3 B = arg^^^r) + i^T^ sm(^^n)} (17) 

When the initial state is a maximal entangled state, i.e., r = 1, the Schmidt coefficients 
are degenerate and the reduced density matrix of the subsystem at t = is 

p a(b \o)= 1 -i. 

Because the system subjects to the bilocal unitary evolution R(t), the reduced density 
matrix of the subsystem at any time t is 

p MB)(t) = Tr B{A) [R A (t) ® R B (t)p(0)R A ^t) ® R B ^(t)] 

(18) 

which means 

I3 A{B) = 0. (19) 

From Eqs. (17) and (19), we can obtain the relation between the geometric phase of the 
subsystems and the entanglement degree. Again the geometric phase of the subsystems 
j3 A ^ is also only a function of the entanglement degree r which can be seen from Eq. 
(17) and is shown in Fig. 2. It can be seen form Fig. 2 that j3 A ^ is a monotonic 
decreasing function of r. If we could measure the geometric phase of the subsystems 
(3 A ( B \ we can infer the entanglement degree r of the protected entangled state. 

4. Conclusion 

In this paper, we have considered two TL ions in a heavily damped cavity. Using 
engineering reservoir scheme, we have obtained a time-independent master equation, 
and then have found a DFS for this master equation, which can be nonadiabatically 
controlled by the system-reservoir parameters. We have achieved a class of decoherence- 
free cyclic evolution of the entangled state. Finally, we have calculated the geometric 
phases of the whole system and the subsystems under the nonadiabatic coherent 
evolution, and have found that there is one-to-one correspondence between the geometric 
phase of the whole system, the geometric phase of subsystems and the entanglement 
degree r . 
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CAPTIONS 

Fig.l The fidelity F as a function of time. 

Fig. 2 The geometric phase of the subsystems /3 A ^ as a function of entanglement 
degree r. 
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